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^ ! Abstract 

a^ : 

• We compute the photon emission rate of an equihbrated, hot QCD plasma 

^ . at zero chemical potential, to leading order in both qem and the QCD cou- 

I pling gs (T) . This requires inclusion of near-collinear bremsstrahlung and in- 

I elastic pair annihilation contributions, and correct incorporation of Landau- 

I ' Pomeranchuk-Migdal suppression effects for these processes. Analogous re- 

I ■ suits for a QED plasma are also included. 

o ■ 
^: 

Ph! I. INTRODUCTION 

^ ' How brightly does the quark-gluon plasma glow? With the commissioning of RHIC and 

^ ■ the future commissioning of the LHC in heavy ion mode, this is becoming an experimental 
• 1— I ! question. It seems appropriate, then, to complete a theoretical analysis of photo-emission 
^ I from a central heavy ion collision at RHIC or LHC energies. 

I The goal of this paper is more modest. We will be concerned only with a fully equilibrated 

quark-gluon plasma, at a temperature high enough that perturbation theory is applicable.^ 
In this context we will compute the photo-emission rate to leading order in the strong 
coupling constant Qs- By this, we mean the spontaneous emission rate for photons of a given 
momentum k, which we will denote as 



z.e(k) = (27rf (1.1) 

Our normalization is chosen so that P^ = / i^e(k) d^h/ (27r)^ is the total spontaneous photon 
emission rate per unit volume. 



^ Since we take the plasma to be thermalized, the spatial size R of the plasma must exceed the mean free 
path for large angle scattering of a quark or gluon; this requires RT ^ 1 at temperature T. 

^ If /^(k) denotes the phase space number density of photons, and = J /-,(k) d'^k/(27r)'^ is the total 



The evaluation of the photo-emission rate has quite a long history |T|-§[ , and was thought 
to have been settled ten years ago when Kapusta et al. ^ and Baier et al. |jTO| independently 



computed the diagrams shown in Fig. |I|, and obtained the same answer. The result of 
Refs. written in a slightly more general fashion than given in those papers, is 



i/e(k) = A{k) [ln(T/moo) + \H2k/T) + C2^2{k/T)] , (1.2) 
where k = |k| is the energy of the photon. The leading-log coefficient A{k) is given by 



A{k) = 2 OEM 1. 



2 1 T^oo 



k 



nf{k), 



(1.3) 



where nf{k) = [exp(A;/T) + l]^^ is the Fermi distribution function, and dp is the dimension of 
the quark representation (so dp = = 3 for QCD). The charge assignment for each quark 
species s is denoted g^, and equals 2/3 for up type and —1/3 for down type quarks. The 
mass moo appearing in Eqs. (|1.2|) and (|1.3|) is the leading-order asymptotic (large momentum) 
thermal quark mass,0 given by [|r^ 



:i.4) 



with Cf the quadratic Casimir of the quark representation. For the fundamental represen- 
tation of SU(iVc), Cf = (Ar2-l)/(2Ar,); for QCD specifically, Cp = 4/3 and = g^T^S. 
Hence, for two-fiavor QCD, the leading-log coefficient is explicitly 



Aik) 



40 7rT2 nf{k) 
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k 



1.5) 



The strong coupling gs (and = S's/^tt) should always be understood as defined at a scale 
of order of the temperature.^ 



number of photons per unit volume, then the evolution of the photon phase space density is given by 
df^{k.)/dt = t'o(k) [1 + ff(k)] — i^a(k) /^(k), where i'c(k) is the spontaneous emission rate and i'a.(k) is the 

corresponding absorption rate. As always, detailed balance relates these rates, i/a(k) = t'o(k) e^^^^"^ . Photon 
absorption can be ignored if the plasma is optically thin; this requires that i?fa 1- For a quark-giuon 
plasma at temperature T, this condition (as we will see) becomes as OfEM RT <C! 1. This should hold for all 
realistic heavy ion scenarios. In contrast, photon absorption can never be ignored in a thermalized relativistic 
QED plasma, since the photon absorption rate is comparable to the electron large-angle scattering rate. 



^ The momentum-dependent thermal mass m{p) is defined by the relation E(jj) — y^p^ + m{p)^, where 
E{p) is the energy of a "dressed" quark moving through the plasma with spatial momentum p (and whose 
zero-temperature mass is negligible compared to T). For a typical momentum of order T, one may expand 
the square root and see that the thermal correction to the quark dispersion relation is proportional to the 
thermal mass squared, AE(ji) = E(j)) — p ~ m{p)'^ / {2p). For p — 0{T) or larger, the momentum-dependent 
thermal mass equals the the asymptotic mass moo, up to irrelevant sub- leading corrections suppressed by 
additional powers of gs- This asymptotic value of m{p) is related to the thermal mass mq, conventionally 
defined as the energy of a dressed quark which is at rest in the plasma, by a square root of two, — Smg. 

''For QED, Cp — 1 and — Hence A{k) equals 2^0^-^ nf{k)/k for an electromagnetic plasma 

containing only electrons and positrons. 



2 




+ 



FIG. 1. Two-to-two particle processes which generate the leading logarithmic contribution to the 
photo-production rate, and were originally believed to give the complete leading order contribution. Time 
may be viewed as running from left to right. 






+ 



FIG. 2. Inelastic processes which, because of near-coUinear singularities, contribute at the same order 
as the two-to-two particle processes. The diagrams on the left represent bremsstrahlung, and those on the 
right are inelastic pair annihilation. Again, time may be viewed as running from left to right. 



The final term C2^2{k/T) in expression (|1.2|) is a function oi k/T witli a finite limit as 
k/T —>■ oo. The papers [P, p!0|] compute C2^2{k/T) only in this limit, and find 



lim 

k/T-^c 



C2^2{k/T) = - 



1 



7e 
2 



+ 



21n(2) , C(2^ 
2C(2) 



+ 



-0.3614902 



(1.6) 



We have also verified this result. 

However, as recently demonstrated by Aurenche et al. |T^-p!^ , the result ( |1.2| ) is incom- 
plete. The bremsstrahlung and inelastic pair annihilation processes shown in Fig. |^ contain 
collinear enhancements which cause them to contribute at the same parametric order in 
coupling, O(aEMtts), as the two-to-two processes of Fig. ^ even for large photon energies 
k^ T. But the calculation of these processes presented in is also incomplete,[]as it fails 
to incorporate an 0(1) suppression of these processes due to multiple scattering during the 
photon emission process, which limits the coherence length of the emitted radiation. This 

The presence of LPM 



is known as the Landau-Pomeranchuk-Migdal (LPM) effect |[T6|-|T8 



which discussed the 



suppression was pointed out in later work by Aurenche et al. |[13|JT^ 
physics involved but did not attempt to make a complete calculation. 

This paper presents a full calculation of the photon emission rate at leading order in g^. 
We evaluate the non-logarithmic two-to-two contribution C2^2{k/T) for general k/T. We 
also compute the rate of photo-production by bremsstrahlung and inelastic pair annihilation, 
fully including the LPM effect. This requires solving a non-trivial integral equation to 
determine this rate. We derived this integral equation in Ref. [jl9[; it is related to one 
written down over 45 years ago by Migdal |T^JT^ for studying electromagnetic energy loss 



'It is also wrong by a factor of 4 p 



3 



of a relativistic particle traversing matter. Related work in the QCD context may be found 
in Refs. [i0-25].| Other discussions of the LPM effect for electromagnetic interactions in 
ordinary matter include Refs. p7|-p9|| . The current paper is a sequel to our earlier work [|19 



but it should be fairly self-contained. 

The outline of the remainder of this paper is as follows. In Sec. ^ we briefly review our 



previous work [jl9l showing why the processes considered here (and no others) contribute to 
the leading order emission rate, and deriving an integral equation whose solution determines 
the bremsstrahlung and inelastic pair annihilation contributions to the emission rate, fully 
including the LPM effect. Sec. |T| discusses, qualitatively, the relative importance of the 
LPM effect in different kinematic regimes. The numerical solution of the integral equation 
determining the LPM contributions to the emission rate is presented in Sec. This section 
is the meat of the paper. Sec. |V| contains the evaluation of C2^2{k/T) for arbitrary photon 
energy k ^ g^T. (If this condition is violated, then thermal corrections to the dispersion 
relations for incoming or outgoing particles can no longer be neglected in these 2^2 
processes.) Our results are combined in Sec. 0, which is followed by a brief conclusion. 

For the convenience of readers interested in just the bottom line, we summarize our 
results here. The complete leading-order photon emission rate may be written as 



with 



i.e(k) = A{k) In {T/m^) + C,,,{k/T) 



Ctot{k/T) = \ In {2k/T) + C2^2{k/T) + Cbrem(fc/T) + Ca„nih(fc/T) , 



:i.7) 
:i.8) 



and the coefficient A{k) given in Eq. ( |L3|) . The function C^otik/T) is the total "constant 
under the log"; it is a non-trivial function of k/T but it is independent of the strong cou- 
pling Qs, whereas ln(T/moo) ~ \n{l/gs) since moo = 0{gsT). Corrections to Eq. ( |1.7|) are 
suppressed by one or more powers of gs- 

The functions C2^2(yk/T), ChTem{k/T), and Canmhik/T) all involve multidimensional in- 
tegrals (or integral equations) which cannot be evaluated analytically but may be computed 
numerically. Individual plots of these functions, as well as their sum, are shown in sec- 
tion 0. Our numerical results for QCD plasmas are reproduced quite accurately by the 
approximate, phenomenological fits 



^2^2 (a;) ^ 0.041a;- 



and 



-0.3615 + 1.01 e~i-^^^, 
0.548 log(12.28 + 1/x) 



0.133 a; 



X' 



3/2 



l + x/16.27 



(1.9) 

;i.io) 



where A^f is the number of quark flavors. The approximation ( |1.9D to the function 6*2^2(2;) 
(which crosses zero) is accurate to within an absolute error of 0.02 over the range 0.2 < 



^References and p^, which predate our work [f9|, derive an integral equation similar to our result. 
The main difference is mat we treat the actual distrioution of moving colored charges in the plasma, with 
dynamical screening, rather than employing a model of static colored scattering centers as in [ pO|pl ] . Rcf. |l9| 
also includes plasma induced corrections to quasiparticle dispersion relations (which do influenceme leading- 
order result), and contains a fairly elaborate diagrammatic power counting analysis to convincingly identify 
all leading order effects. 
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k/T < 50, but this form fails to be a good approximation at smaller fc/T; 6*2^2 (a;) has 
a finite small x limit, rather than growing as 1/x. As discussed in section pT|, the form 
( |1.10| ) for the near-coUinear contributions builds in the right parametric large and small k 



asymptotic behavior, up to logs, and is accurate to within a 3% relative error over the range 
0.2 < fc/T < 50 for A^^f from 2 to 6.Q 

Consider, for the sake of discussion, a value of ctg = 0.2 which corresponds to ln(T/ moo) = 
0.09. In this case, it turns out that the old results of Baier et al. and Kapusta et al. [Q, 
which neglect near-coUinear bremsstrahlung and inelastic pair annihilation, are within a 
factor of two of the correct leading-order rate for photon momenta in the range 2.5 < k/T < 
10. Bremsstrahlung becomes the most important process for energies below this range, 
and inelastic pair annihilation makes a large relative contribution to the rate above this 
range (where the total rate is small because of exponentially falling population functions). 
The LPM effect suppresses both bremsstrahlung and pair annihilation processes, but the 
suppression is not severe (35% or less) except for soft bremsstrahlung with k < 2T, or 
very hard pair annihilation with k ^ lOT, where the LPM suppression is significant. For 
parametrically small photon energy A; < T (but with A; > eT and A; > ^^^Tln^-^), the 
LPM effect changes the parametric behavior of the emission rate. One finds, 

dV^ock-^/^dk, (Lll) 

[ignoring log(T/A;) factors], rather than k~^ dk which is the result if one ignores the LPM ef- 



fect. Therefore some results in the literature |12,26] concerning soft, on-shell photo-emission 



from the quark gluon plasma, ignoring the LPM effect, are wrong parametrically. Similarly, 
for extremely hard photon emission, the rate behaves as fc^^^ e~''^^dk, rather than fc^ e~'^/'^dk 
as one would obtain by neglecting the LPM effect. For more detailed results see the figures 
in Sec. |VI[ 



II. BREMSSTRAHLUNG AND INELASTIC PAIR ANNIHILATION 

The photo-emission rate is related to the imaginary part of the electromagnetic current- 
current correlator. We show in that a correct treatment of the near-collinear 



bremsstrahlung and inelastic pair annihilation processes which contribute to the leading- 
order emission rate requires a summation of all ladder graphs of the form shown in Fig. |^. 
(Self-energy insertions on all propagators, including the leading order imaginary parts, are 
implicitly resummed.) As indicated in Fig. ^ these diagrams can be interpreted as gen- 
erating various contributions to the squares of scattering amplitudes. It is necessary to 
sum all ladders because there can be interference between emission at one scattering and at 
subsequent scattering events. This is the basic origin of the LPM effect. 



^ It should be emphasized that throughout our analysis Jl9] we assume that the photon wavelength is much 
smaller than the large angle mean free path for quarks. TMs requires that the photon energy be parametri- 
cally large compared to g^T\n.g~^. We also req uire k 3> m^, where is the asymptotic (electromagnetic) 
thermal mass of the photon, given in Eq. (pTsf). This latter condition ensures that plasma corrections to 
the photon dispersion relation are small and do not significantly reduce the photon velocity below c. Given 
the physical value of q;em, this implies that our results for soft bremsstrahlung can only be trusted down to 
about k K, 0.2 r (for QED with only electrons and positrons) or fc w 0.3 T (for three flavor QCD). 
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FIG. 3. Typical ladder diagram contributing to the electromagnetic current-current correlator. All 
such ladder diagrams must be summed to determine the leading-order bremsstrahlung and inelastic pair 
annihilation rates. Resummation of self-energy insertions on all propagators is implied. 




FIG. 4. Cut ladder diagram with N gluon "rungs" may be interpreted as the interference between 
amplitudes for photon emission before and after N scattering events. 



The total photon emission rate is given by an integral over all particle momenta of the 
emission rate for that particle momentum. Factoring out the same coefficient A{k) which ap- 
pears in the leading logarithmic result ( |1.2|) , we find that the contribution of bremsstrahlung 
and inelastic pair annihilation processes to the leading-order emission rate is 



'(k) = A{k) / dp 



pf + {P\\+kf 



mi 



nf{k) 



(2vr)^ 



Re {2pj 



f(P±,P||,fc)}- 
(2.1) 



Here p\\ and are the components of the momentum p parallel or perpendicular to the 
photon direction, respectively. The contribution from the integration region —k < p\\ < 0, 
in which py and {p\\+k) have opposite sign and l—nf{p\\) = nf{\p\\\), should be interpreted 
as inelastic pair annihilation, while the contributions from > or p\\ < —k represent 
bremsstrahlung off a particle of energy p\\+k (when p\\ > 0) or an anti-particle of energy 
-P\\ (for p\\ < -k). 

The final term Re{2p^ ■ f (pj_,p||. A;)} is proportional to the amplitude squared for a 
particle of momentum p-|-k to emit a photon of momentum k, and is the non-trivial part 
of the calculation. The function f(px,p||, k) is the solution to an integral equation, derived 
in [O, which accounts for the partial coherence between multiple scattering events. This 



integral equation is 



2p^ = z5E(p^)f(p^)-f 



(27r)2 



C(q±) f (px) - f (p±+q±) 



(2.2) 



where the dependence of f (pj_,p||, /c) (as well as 5E) on p\\ and k has been suppressed for 
notational convenience. The collision kernel C(qx) is given by 



6 




212 



(2.3) 



where g = |q| = [q5_ + Q'li^]"'^''^, tti-d is the (leading order) Debye mass, and Iii^{q^,q) and 
nT(9°, Q') are the standard hard thermal loop longitudinal and transverse gluon self-energies 
Pm , respectively; their explicit form is given in Section In the first term on the right-hand 
side of Eq. ( p.2| ), 5E denotes the energy difference, on shell, between the two relevant states 
of the system, in the limit where p\\ is large compared to |p_l| and moo- For bremsstrahlung, 
this is the energy difference between a quark of momentum p plus a photon of momentum 
k, and a quark of momentum p-|-k, which explicitly is 

Pi + "41 + ^ ■ (2.4) 



For pair annihilation, 5E is the difference between a quark plus anti-quark with momenta 
— p and p-|-k, and a photon of momentum k; this gives exactly the same result (|2.4|) given 
our convention about the sign of p\\. This energy difference generates a phase difference 
between the amplitudes in Fig. ^, which is manifested as an imaginary term in Eq. ( ^.21) . 
This phase difference will eliminate coherence in the emitted radiation beyond a time scale 
of order 1/SE, which is often referred to as the formation time of the photon. The 
contribution to 6E accounts for plasma induced dispersion in the photon propagation; is 
the asymptotic thermal photon mass, whose square is half the electromagnetic Debye mass 
squared of the plasma. 



d^Y^^ne'T^. (2.5) 



Since we assume oem ^ «s, we may neglect this electromagnetic dispersion correction 
when considering photo-emission from a quark-gluon plasma; but when we consider the 
emissivity of an electromagnetic plasma, which we will also compute, retaining this term 
will be essential. 

Note that the collision kernel (|2.3| ) is rotationally invariant (in the transverse plane). 
Consequently, the solution to the integral equation ( |2.2| ) must equal p_L times some scalar 
function of p^, 

f(p±) = P±x(pl)- (2.6) 



Hence, although written as a vector equation, the integral equation ( p.2|) is trivially reducible 
to a scalar integral equation. For later convenience, also note that this integral equation 
may be written more abstractly as 2pj_ = [i 6E + C ] f , if one defines C as a linear operator 
acting on functions of p_L, whose explicit form is given by the second term in Eq. ( p. 21) . 
Further discussion of equation ( |2.2| ), and its derivation, may be found in Ref. p9[ . 
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III. LPM EFFECT: ROUGH ESTIMATES 



Before solving equation (|2.2|) numerically, it is useful to examine in what kinematic 
regimes the LPM effect will be important and to estimate how large its effect can be. 
Neglecting the LPM effect altogether means ignoring any interference between successive 
scattering processes. This corresponds to solving Eq. ( p.2p under the assumption that 5E 
is much greater than the collision term which follows it.| At leading order in this "approxi- 
mation," f(p±) is then pure imaginary. 



Imf(p_L) 



2pj 



no LPM 5E{j>i.) ' 



(3.1) 



The next order arises from substituting this result into the collision term, thereby determin- 
ing the leading real part of f in this approximation. 



Ref(p±^ 



no LPM 5E(j>±) J (27r 



C(qi 



2p± _ 2(p^+q^) 



(3.2) 



Inserting this form for Ref(px) into Eq. ( ^.11 ) reproduces the results of Aurenche et al. 
(except for the overall factor of 4 error made there [l^). Given the form (|2.4|) of 5E, and 
neglecting m^, one may see without performing any calculation that Ref(p_L)|noLPM will be 
ap|j and k independent function times p||^(j9||+A;)^/A;^. Therefore, in this approximation, the 
photon spectrum will look like 



Z/e(k) 



oc Aik) I dpw 

no LPM ^ ' J-k/2 " 



+ {P\\+kf 



nf{k+pii) [l-n/(p||)] 
nf{k) 



(3.3) 



with a ^-independent constant of proportionality. For soft photons {k <^ T) the dominant 
contribution comes from ~ T, and the result, ignoring the LPM effect, grows as T^/Zc^ 
relative to the emission rate due to only 2^2 processes. At large k ^ T the dominant 
contribution comes from the whole range —k/2 < p\\ < over which the population function 
terms are approximately unity. This yields a result (again, ignoring the LPM effect) which 
grows as k/T relative to the 2^2 emission rate. Therefore, ignoring the LPM effect, 
we find that bremsstrahlung dominates the emission rate for k <^ T, while inelastic pair 
annihilation dominates the rate for k^T. 

We must now consider whether, for any range of k, it is actually permissible to neglect 
the LPM effect. Using the form ( p.4| ) of 5E, one may easily see that if kT/\p\\ {p\\+k)\ is 
sufficiently large then \5E\ will also be large for all p^ and an expansion based on \5E\ ^ C 
will be valid for all p^. However, in this regime, Ref is always small, and so is the photo- 
emission rate. In both cases where we have just found a large photo-emission rate (ignoring 
LPM suppression), the quantity kT/[p\\ {p\\+k)] is small, and hence the inequality 5E ^ C 



^Note that 5E does not vanish for any value of p_L and p\\ provided the inequahty 2mao > is satisfied. 
This holds in QCD because of our assumption that e is much smaller than gg. It is also satisfied for electron 
and electron plus muon QED plasmas. If this inequality is violated, then a photon propagating through 
the plasma could dec ay into a fermion-antifcrmion pair without any associated inelastic scattering, and the 
resulting pole in ( |3.l[) would produce a real contribution to f (p±) proportional to pj^ S{SE{p±)). 
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only holds for large p_L. Therefore, the estimates we have just made are too optimistic, and 
there will be substantial LPM suppression in exactly those regimes where the inelastic rate 
is dominant. To estimate how large the suppression is, we may solve for f(p_L) at large p_|_, 
where the large 5E approximation is valid, and determine where the approximation breaks 
down. At large p_L, we have 



Imf(pj 



2pj 



I asymptotic 



5^(P±) 



and 



2p\\ {p\\+k) 



2rl 

Pi 



(3.4) 



Ref(p±: 



I asymptotic 



5^(P±) 



n TvYi f I 

±111 1 I asymptotic 



"2p|| {j)\\+k)' 


2 


k 





(2vr)2 



C(q±) 



2p± 2(p^+q^) 



Pi 



(3.5) 



This integral may be evaluated in closed form. Since the kernel C(q 
ant, the integral in Eq. (|3.5|) may be rewritten as 



rfjqj 
27r 



|q±|C(|q±| 



2p± 

Pi 



2pj 

Pi 



(p±+q±)2 
l) is rotationally invari- 

(3.6) 



where {g{p±))q^ denotes the average of a function g over the circle of radius q± = |qj_| 
centered at p_L. Next, note that (both components of) p±/pi are harmonic functions 
(away from p^ = 0), and recall that the average of a harmonic function around a circle 
equals the value of the function at the center of the circle, provided no singularities are 
enclosed. Consequently, the integrand of Eq. ( p. 61 ) vanishes for |q_|_| < |p_l|. For |q^| > |p_l|, 
elementary complex analysis shows that (p±/pi)„ vanishes.^ Therefore, 



Ref(p_L)|^3y^pt„^i^ 



1 

7r 



2p|| {p\\+k) 



k 



P± 



\p±r ^ipxI 



C(|q±|)|q±|rf|qj 



(3.7) 



Since pj_ has been assumed to be asymptotically large, and |qj_| > |p±|, one may completely 
neglect the self-energy corrections in the explicit form (|2.3| ) of C(qx). The resulting integral 
is then straightforward, leading to 



Ref(p±) 



I asymptotic 



2p\\{p\\+k) 



k 



ttT 



2pi 
|P±I 



(3.8) 



These large p^ asymptotic forms are valid when 6E ^ C, or equivalently when 
Re f(p_L) [asymptotic ^ f (p_L ) | asymptotic- The cxpausiou wlll havc broken down when the 
real and imaginary parts become comparable. This occurs for 



^If one relabels p_L as r, and q± as R, then the reasoning just given is easily seen to be just a two dimensional 
version of Newton's classic result that the gravitational force due to a uniform spherical shell of radius R 
vanishes inside the shell, but outside the shell is the same as the force due to a point mass at the center of 
the shell. 
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p1 ~ 



Pj 



crossover 



lp\\ {P\\+k) 
kT 



(3.9) 



Below this crossover point, Ref(p_L) will be smaller than the asymptotic estimate. The 
region ~ P±| crossover will, in fact, provide the dominant contribution to the final cPp± 
integral in the rate (|2.1| ). Hence the size of this integral [ignoring factors of ln(A;/T)] may 
be estimated as 



/; 



'^(pI) R-S f I asymptotic ' P± 



I crossover 



P\\ 



k 

P\\iP\\+ k) 
kT 



crossover 



-1 3/2 



mo moo T . 



(3.10) 



Using this result, we find that, for bremsstrahlung, the small k emission rate grows, 
relative to the 2 2 rate, only as (T/fc)^/^, not (T/ky as was found ignoring the LPM 
effect. Similarly, for large k the inelastic pair annihilation rate grows, relative to the 2^2 
rate, only as {k/Ty/"^, not {k/T) as was found when ignoring LPM. The dependence on 
TTT-D is also simple; relative to the 2^2 rate, the inelastic contribution grows linearly in 
the Debye mass m-Q. Therefore, inelastic processes become relatively more important as we 
increase the number of species in the plasma (which increases ttid), but only as the square 
root of the number of species. 



IV. SOLVING THE INTEGRAL EQUATION 

In order to solve the integral equation (|2.2| ) which determines f(p±), we will use an 
adaptation of the variational approach we previously exploited in pT| , |5^ . But first, it will 
be convenient to scale out the dimensionful parameters from Eq. ( p.2|) . To do so, we define 
dimensionless transverse momenta, 

p± = — , q = — , (4.1) 

mD mo 



5E(px) ^ — 5E{p^) , f (p^) = f (p^) . (4.2) 
mfv 1 



and correspondingly define 

T_^r^/__ N r/~ N _ mD 

It will also be convenient to define the dimensionless ratio of thermal masses 

_ <o 3Cf [ 3/(4iVf) , QED; 

ml 4{TA + NiTF) \2/(6 + iVf), QCD . ^^""^ 

Here Ta and Tp are the trace normalizations of the adjoint and quark representations, 
respectively. For QCD, Ta = iVc = 3 and Tp = 1/2, while for QED, Ta = 0, and Tp = 1. 



Substituting these dimensionless variables into Eq. (^TTP , and noting that the argument 



of the p|| integral is invariant under p|| —k — p\\, that equation becomes 
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^LPM^j^) ^ Ak)T 1^ ^^^^ fp||2 + [pw+kf] nf{k+p\0 [l-nf{p\0] 



K J-k/2 



nf{k) 



(2p±, Ref 



(4.4) 



where we have introduced a natural inner product on the space of (vector valued) functions 
of transverse momentum, 



(fi , f.) 



(2vr)2 



fl(P±) ■ f2(p±) • 



(4.5) 



Once again, the contribution from < is inelastic pair annihilation, while the contribution 
from p|| > is bremsstrahlung. The dimensionless form of the integral equation for f is now 



where 



and 



with 



5E = { 



2p± = 


'i6E + c] f , 


r kT 






2p|| {k+p\\) 




kT 




T 

^ 4k 


. 2p\\ {k+p\\) 



QCD; 
QED, 



(Cf)(pj 



{2ny 



m±) f (p±) - f (p±+q±: 



(4.6) 



(4.7) 



(4. 



C(qi 



dq° 5{q\\-q 



[1 - (g7g)2i2 



. z 2 + - 

1 \ q^- nL(gO, q) (g0)2 _ ^2 _ n^^^o^ ^) 



(4.9) 



Here gy = gy/mD, q'^ = q^/mD, and q = y^ql + g|- The explicit forms of the thermal gauge 
field self energies are 



nL(g°,g) 
nT(g°,g) 



-1 + 



2q 



In 



ITT 



z0\2 



2q 



nf2 



+ 



e - iff 



4g3 



In 



ITT 



(4.10) 
(4.11) 



These forms are valid for q > \q^\, which is enforced by the delta function in ( |4.9| ) which 
identifies with gy. Both 6E and C should be viewed as linear operators on the vector 
space of functions of transverse momentum. Both operators are Hermitian — in fact, real 
and symmetric — under the inner product ( ^.5| ). 
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A. Variational approach 



We are unable to find an exact solution of Eq. ( [4 .61) for arbitrary values of k and 
kT / p\\{k+p\\) . Instead we must use some approximative, but accurate, method. We want 
a method which can be made arbitrarily accurate as computational effort is increased, and 
which in practice will have relative errors of less than 10~^. This goal can be met with an 



adaptation of the variational approach used in our previous work [pT| , |32 
First, separate the complex equation 



into a coupled pair of real equations. 







C 


-se' 




"Ref ■ 







_-5E 


-C 




Im f 



(4.12) 



The linear operator represented by the 2x2 matrix is real, symmetric, and non-singular. (In 
fact, eigenvalues of this operator are bounded away from zero by at least \k kT / [2p\\{k+pi\)]\.) 
Consider the functional 



2p±,0 



Ref 
Imf 



Re f , Im f 



C_ -6E 
-6E -C 



Ref 
Im f 



(4.13) 



It possesses a unique extremum, at which Eq. 
extremum is 



is satisfied. Moreover, its value at the 



Q. 



extremum 



i(2p^, Ref 



(4.14) 



which is precisely half the quantity appearing in the final integral ( ^.4| ) for the emission 
rate. However, the extremum is a saddle-point, rather than a maximum or minimum. 
Nevertheless, if one uses a sufficiently flexible Ansatz for f , then the extremum under that 
Ansatz will be very close to the true extremum and the determined value of (2p^, Ref) 
should be very accurate. But because the extremum is a saddle-point, a variational answer 
based on a finite basis set does not automatically yield an upper or lower bound on the true 
answer, nor is there any guarantee of monotonicity as the size of the basis set is increased. 

As noted earlier [c.f. Eq. ( |2.6| )], rotation invariance implies that f(p±) equals p^ times 
some scalar function x{P'±)- So, as an Ansatz for f, we will choose a finite linear combination 
of terms of this form. 



Ref(pi^) = p_ 



m=l 



Imf(p^) = p_ 



m=l 



Am) I ~2 



(4.15) 



There is no need to choose the same basis of trial functions for the real and imaginary parts. 

it is in fact sensible to 



Given the differing asymptotic behaviors [c./., Eqs. ( |3.4| ) and 
use different basis sets for the real and imaginary parts. There is also no need to require 
that the number of trial functions for the real and imaginary parts, iVr and A^;, be the same; 
but in practice we will do so for simplicity. Given this Ansatz, the functional Q[f\ becomes 
a quadratic function of the coefficients {am} and {6m}, 



Q[f] = Q{{am], {hm}) = ^a'mSra-]:^\^ 



-6Enm ~Cmn 



(4.16) 



12 



where 



(4.17) 



The extremal value of Eq. 16|) is given in terms of these components by elementary linear 
algebra as, 



Q 



extremum 



where 



S 







M 



-SE 



(4.18) 



(4.19) 



are the indicated A^r + component vector and (symmetric) matrix. 

It remains to choose the trial functions and perform the integrals. The choice of trial 
functions should be informed by examining the general behavior of the integral equation 
]6|) we are trying to solve. At small p]_, x(p1) i^^iy be approximated as a polynomial in 



p\. At large p\, we have already determined its behavior; from Eq. ( p.4|) and Eq. (|3.8| ) one 
sees that the imaginary part of x(Pj_) must go to zero as l/pj_ while the real part falls as 
l/{p\)^. We may choose trial functions so that every trial function decreases at large 
at least as fast as required; at least one trial function must fall with the correct asymptotic 
power. A nice set of trial functions which have the appropriate form is 



(pIM)--i 



;i+piM) 



m = 1, ...,N, , 
m = l,...,Ni. 



(4.20) 
(4.21) 



We have not specified the constant A, which should be chosen to be close to the value of p']_ 
where f is varying in the most significant way. For instance, it can be chosen to be the value 
of p\ where Re f reaches half of its asymptotic value, or the value at which the integrand 
in (2p±, is peaked. The only difficulty is that the value of p^ where f shows significant 
behavior depends on k, k, and p\\. We are interested in different values of n and k, and p\\ 
is integrated over in Eq. ( [4.4|) . Therefore a single value of A will not give a good set of trial 
functions for all values of k, k, and p\\. Our solution is to determine all the matrix elements 

{Cmn} and {6Emn} for a variety of values of A. Then we determine which value of A to 
use, for a particular set of parameters [K,,k,p\\], based on a criterion which optimizes the 
sensitivity of the trial functions to the most important range of p\. Various specific criteria 
could be used; we chose the following. We want the integral 



(2p^,f)=— / Rex(pi)pid(pi 
^ ^ /vr JO 



(4.22) 
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to be accurately determined. Map the integral to a finite interval by defining 



y ^ Ain ZV I ' (4.23) 
A.^' + |p±| 

and rewrite the integral in terms of 

(2px.f)=i|^X-(V/(l-=/)V!'- ("■24) 

For each value of A, we solve for the variational coefficients {a™,, &m} and determine the 
resulting function x{P\)- We then choose the value of A for which the maximum of the 
integrand of Eq. ( [4.24| ) occurs closest io y = 1/2, which is where our basis functions have 
the most sensitivity.|^ This ensures that we are always using a set of trial functions which 
are suitably flexible in the region most important for evaluating the integral that we care 
about. 

The final integral over p\\ in Eq. ([4.4|) is one dimensional, and performing it by numer- 
ical quadrature is elementary. Doing so completes the calculation of the contribution of 
bremsstrahlung and inelastic pair annihilation to the photo-emission rate. 

B. Integrals for Cmn 

The integrals determining 6Emn and Sm are elementary for our choice of trial functions 
and can be done in closed form. One finds 

Sm = 2A^ K{m,N,) , (4.25) 

5Emn = ^ , , [AK{m+n,N,+N,) + r]K{m+n-l,N,+Ni)] , (4.26) 
2p||(/c+p||) 

where K{j,N) = j\ (A^-j)!/[47r (A^+1)!] and r/ = k for QCD, or k + ^ p\\{k+p\\) / k"^ for QED. 
The collision integrals {Cmn} cannot be evaluated analytically, but may be reduced to a 
numerically tractable form as follows. The integral we need is 



^mn--,l (2^)2 7 (27r)2^('^^^^ 



Px0^™^(p1) - (p±+q^)0^"^^(|p±+qj^)j ■ [p±0("Hp1) - (p±+q±)0('^)(|p^+qj^^ 

(4.27) 

Here we have used the symmetry of the collision operator to write the integrand in a form 
manifestly symmetric under m <-> n, which also makes apparent that the collision operator 



"'The maximizing value of y sometimes depends on A such that it crosses y — 1/2 not only at A 



y'l + {k[k+p]/kT), but also at a very large value of A. This second solution is spurious and should be 
discarded. 
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is positive definite on all normalizable functions f . One overall angular integral may be 
performed, and the remaining integrals become 



Cn 



327r2 



C{q±)Ip{q±) , 



J-TT In L 



+ |p±+q± 



|2^M|'li=, 



(|p±+q±l 



(4.28) 



- (Pl+P±q± cose) [0("nip±+q±n <P^'"\pI) + (m ^ n)] | , (4.29) 

where in the last equation (p±+q±)^ is shorthand for p\ + q± + 2p±q± cos 6. Although we are 
not able to evaluate either Ip{q±) or C{q±) analytically, since no more than three integrals 
ever nest it is not difficult to perform all the integrals by numerical quadrature, using 
adaptive mesh refinement, and achieve better than a part per million relative accuracy.0 

We have evaluated the quantities {C^n} {^mn}' several values of A in the range 
0.6 to 20, using = N,. = N ranging from 2 up to 16. Reducing from 10-16 to 6-10 
changes our results for the integral (|4.4| ) by less than 50 parts per million for all values of 
k and k which we consider [k > T/10). Using only 4-6 trial functions causes errors of less 
than 0.2%, while using 2-4 trial functions allows errors of a few percent. We conclude that 
this variational approach, with a modest set of test functions and fiexible selection of the 
overall scale A, yields a very accurate determination of (^p±,i), and likewise for the final 
bremsstrahlung and inelastic pair annihilation rates. 



V. 2 ^ 2 PROCESSES BEYOND LEADING-LOG 



As discussed in the Introduction, bremsstrahlung and inelastic pair annihilation processes 
are not the only processes which lead to photo-emission at order agaEM- There are also the 
2^2 particle processes depicted in Fig. 0. Neglecting for a moment the subtleties associated 
with the kinematic region where the internal propagator becomes soft, and ignoring the 
(small) corrections to the external state dispersion relations due to both thermal effects and 
explicit mass terms, the contributions to the photon production rate from (gluon-photon) 
Compton scattering and quark- ant iquark annihilation are 



pCompton 

7 



T Annihilation 



d^p' ^^^^^ 5\p,+p'-k,-k'^) 16e\Es <il 



(27r)i2 IQpkp' k' 

X dp Cp 



X 



nf{p) Tibip') [l-nf{k')] 



■ 7 



i2nyn6pkp'k' (2-) ^ (P.+P.-k.-K) 16e (E. 
X dp Cf gl 



X 



(5.1) 



(5.2) 



^^When trying to achieve high accuracy, there can be some difficulty in the 9 integral near — n when 
p_L is very close to q± and both are large. This can be eliminated by noting that the integral has a factor 
of two redundancy; the integral we want is twice the integral where 9 is restricted to values for which 

(p±+q_L)^ >pI- 
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respectively, with p and p' the incoming momenta, k the photon momentum, k' the other 
outgoing momentum, s, t, and u the usual Mandelstam variables, and p = \p\ = etc. 

These expressions must be corrected when the momentum on an internal propagator 
becomes small, because thermal effects then significantly modify the matrix element. This 
occurs for the u/t and s/t terms when t is of order g^T^. Therefore it is best to arrange 
these integrations so that |p— k| appears as an integration variable, and to separate off the 
soft (small |p— k|) region for special treatment. Such a separation and treatment of the soft 
region was done correctly in Refs. PJTO[|. Outside this soft region, Eq. ( |5.1|) and Eq. ( ^.2|) 
are correct at leading order in g^. The large |p— k| region was treated in Refs. 0,T^ only in 



the limit T, but we will not make any such restriction in our evaluation here. 

Because the population functions take a simple form in the plasma rest frame, we choose 
to work in terms of plasma frame quantities. We first consider the term in Eq. involving 
the t/s matrix element. This term requires no special handling of the soft region, and is 
most easily done in the "s-channel" parameterization of Ref. Defining uj = p+p' and 

q = I p+p' I, the photo-production rate resulting from the t/s term is 

n-^^^^^l^r-r^/;*/;-r<^^ . (5.3, 







X nf{p) nb{uj—p) [l—nf{uj—k)] x 

X |g2 - {2p-Lj){2k-Lj) - cos(/)y^[g2-(cu-2p)2][g2_(,^_2A;)2] | . 



The integration is trivial and just eliminates the cos0 term. To determine the energy 
dependence of the photon spectrum we want the result at a fixed value of k, which is 

with 



00 



n 



(p) nh{u-p) [l-n/(a;-/ 



hik) ^ — — du , ^ _ dp ^^"^ , ^ q' - {2p-u){2k-u 



(5.5) 



2'K'^T'^Jk J\2k^u\q^J^ kufik) 



We have factored out the same leading logarithmic coefficient A{k) defined in Eq. ( p..3| ), in 
such a way that Ii{k) directly gives a contribution to C2^2{k/T). 

Next we turn to the terms involving s/t and u/t matrix elements. For these terms, let 
UJ = {k—p) and q = |k— p|. Introduce a separation scale q^: satisfying ggT ^ q^ and, as 
mentioned above, separate the soft (g < q^) contributions from the remainder. The emission 
rate from the soft q < q^ part of the s/t and u/t terms requires the inclusion of thermal 
corrections on internal propagators, which modifies the matrix element from what is written 



Our labelinK of external momenta differs from our previous work [pil|32|l. The incoming momenta are now 

3^ the incoming momenta 

were (p,k) and outgoing momenta were (p',k'). The current choice assigns the photon momentum fc, 
consistent with Ref. ||f9| and the previous section. 



labeled (p,p') and the outgoing momenta are (k, k'), whereas in the papers |31 
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in Eqs. and (|5.2|) . However, one can make q <^ T approximations which render the 

integrals tractable. The q < q^, result was found in Ref. [10| and i^ 



dk 



k^A{k) 



27r2 



In 



In 2 



(5.6) 



The emission rate from the remaining q > q^ part of the s/t and u/t terms is 



-p[u/t,s/t] 

7 



dFC^gy {Esq's 



q>q* 



dq 



duo 



(27r)6 Jg, j-q jin^ JO 

X <! nf{p) nb{p) [l-nf{p-uj)] {2p+uj){2p ~uj) + q^ 



# COS( 



+ nf{p)nf{p') [l+ni,{p' -uj)] {2p+uj){2p' -uj) -g^ -#cos0 



+ 



(5.7) 



where we have not written the coefficient on the cos (j) pieces because the dcf) integral simply 
removes these terms. We wish to extract the part of this rate corresponding to fixed k = p+uj. 
The result is 



^YNt,s/t] 



dk 



q>q* 



fcU(fc) 

27r2 



h{k,q^) 



(5.8) 



with 



h{k,q* 



1 



27r2T2 



oo fJ^Q /•min{q,2k—q) poo 

— I duo \ dp' X 

q, q^ J-q 

\nf{k-uj) ribip') [l-«/(p'-^)] 



+ 



knf{k) 

nf{k-uo) nf{p') [l+nb{p' -uo)] 
knf{k) 



{2k-uo){2p' -uo) + q^ 
\2k-uo){2p -uo) - q' 



+ 



(5.9) 



The complete contribution to the emission rate from these 2 
given in Eq. (|1.2|) with 



2 processes has the form 



C2^2{k/T) 



lim 



Ji(fc) + hik, q,) + ln(g,/T) + \ ln(2T//c) - 1 



(5.10) 



For q <^ T, the uo and p' integrals in Eq. ( |5.9|) may be performed analytically and yield 
/2(A;,g*) ~ Iq^dq/q = — Ing* + const., which con&ms that the g* — *• limit in the full 
result ( |5.10|) is well behaved. If one assumes that k ^ T, then the remaining integrals in 
Eqs. ( f).5\j and ( ^^ can also be evaluated in closed form, to leading order in T/k, yielding 
the asymptotic result for C2<~+2 presented earlier in Eq. ( |1.6| ). For general k T, one must 
perform the triple integrals in Eqs. (|5.5| ) and ( ^.9|) by numerical quadrature. To obtain 
accurate results, one should use several values of g^, and perform an extrapolation to small 
g* to remove 0(g*/T) corrections. This we have done. The integrals are not numerically 
demanding. 
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Eq. (17) of Ref. [nOl leaves one un-evaluated integral, which when performed gives this result. Ref. 



computes the {q 



< contribution, which can also be converted to the result quoted. 
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VI. RESULTS 



The complete leading-order spontaneous photon emission rate may be written as 

z/e(k) = {2Txf ^ = A{k) [in (T/moo) + CUk/T) ] , (6.1) 



with the coefficient A{k) given in Eq. ( [L.3|) and 

Ctotik/T) = i In {2k/T) + C2^2{k/T) + C^.^k/T) + C^^^^{k/T) (6.2) 

the total constant under the log. The contribution C2^2{k/T) is determined from the sum 
of Eqs. ( ^.5|) and (p.9|). It is independent of k = ml^/m'^, and is presented in Fig. |, which 
also shows the asymptotic value ( |1.6|) derived in Refs. [P,p!0| . Note that the asymptotic result 
provides quite a good approximation provided k > ST. 




Exact 



Asymptotic 



-0.4 - 



2 4 6 8 10 12 

Photon energy k/T 

FIG. 5. Contribution C2^2{k/T) to the constant under the log arising from 2^2 particle processes 
(gluon-photon Compton scattering and auark-antiquark annihilation), as a function of k/T. The dotted 
line is the asymptotic limit, derived by ]g|,p^. The result is quite close to the asymptotic value already at 
k/T = 3; the visible but tiny difference between the curve and the asymptotic value is due to a 1/fc tail with 
a very small coefficient. 



The functions Cbj.em{k/T) and Cannih(^/^) are k dependent. We show Cbrem(^/^) [de- 
fined as the contribution to the integral (|2.1| ) from the regions p\\ > plus p\\ < —k], and 
Canmh{k/T) [the —k < j9|| < contribution to Eq. (|2.1| )] separately in Fig. |], for a number of 
physically interesting values of k. The values shown were selected to equal the value for an 
Abelian electron plasma (k = 3/4), an electron plus muon plasma {k = 3/8), and a quark- 
gluon plasma with 2 to 6 flavors {k = 2/(6 + A^f)). For the Abelian plasmas, the photon 
dispersion correction is included in 6E [c.f. Eq. ([4.7|)], which suppresses soft bremsstrahlung. 

The effect of LPM suppression on the bremsstrahlung and inelastic pair annihilation 
rates is shown in Fig. 0, which plots the ratio of the actual bremsstrahlung and inelastic 



18 




2 4 6 8 10 0.2 0.4 0.6 0.8 1 



Photon energy k/T Photon energy k/T 

FIG. 6. Left panel: Contributions Chicn\{k/T) and Ca.nn\\\{k/T) to the constant under the log aris- 
ing from near-coUinear bremsstrahlung and pair production. The curves rising at small k/T are plots of 
Cbrcm(fc/T), while the curves rising at large k/T are plots of Cannih(fc/T)- In both cases, the curves from 
bottom to top are for k = 3/4, 3/8, 1/4, 2/9, 1/5, 2/11, and 1/6, corresponding to electron plasma, elec- 
tron plus muon plasma, and QCD with 2, 3, 4, 5, or 6 quark flavors, respectively. In the bottom two 
curves, we include the plasma induced photon dispersion, relevant for QED plasmas, which suppresses the 
soft bremsstrahlung rate. Right panel: Magnified view of the soft photon region of the bremsstrahlung 
contribution. 

pair annihilation rates to the corresponding rates neglecting the LPM effect, for the case 
of iVf = 2 QCD. For most momenta, it is evident that the LPM suppression is a rather 
modest 35% or less. The exception, for the range of photon momenta shown {k < 10 T), 
is bremsstrahlung with k < 2T. Here the LPM effect is important. It also makes a large 
relative suppression to the rate for inelastic pair annihilation at very large k/T, but in 
practice this means k > 20 T, where the statistical factor nf(k) ~ e"'^/"^ ensures that 
virtually no photons are produced. 

The complete constant under the log, Ctot{k/T), equal to the sum of contributions shown 
in Eq. ( |6.2|) , is plotted in Fig. |^. Over a wide range of photon energies, T < k < 6T, one 
sees that Ctot{k/T) lies between about 1 and 2. 

Finally, Fig. ^ shows the full leading-order photon emission rate, together with the 
bremsstrahlung, inelastic pair annihilation and 2^2 contributions, for two-flavor QCD 
with as = 0.2. The left panel shows the production rate dT^/dk, while the right panel 
shows the rate weighted by a factor of the photon energy. One sees that bremsstrahlung 
is the dominant production mechanism for not-so-hard photons with k < 2T. As discussed 
in section |T|, the ratio of the bremsstrahlung rate to the 2^2 rate scales, at small k, 
as {T /k)^/'^ (up to logs), a factor of \Jk/T less strongly than if there were no LPM effect. 
At intermediate values of fc/T > 4, where the most experimentally detectable photons are 
probably produced, the 2^2 processes generate the largest contribution, comparable to 
the sum of the two near-collinear inelastic processes. At large k ^ T, the inelastic pair 
annihilation emission rate grows, in comparison to the 2^2 rate, almost linearly with k/T 
up to quite large k. As shown in section the ratio of these rates ultimately grows only 
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2 4 6 8 10 

Photon energy k/T 



FIG. 7. Relative importance of the LPM effect. The figure shows the ratio of the true bremsstrahlung 
and inelastic pair annihilation rates to the rates which result if one neglects the LPM effect in the treatment, 
for the case of k = 1/4 which corresponds to two flavor QCD. The LPM suppression is rather modest, 30% 
or less, except for bremsstrahlung photons of energy k ^ 2T. These plots show very little change for other 
values of the mass ratio k. 



T , I , , , I , , , I r- 

Bottom to top: 

QED with 1-2 species (dashed), 

QCD with 2-6 flavors. 




2 4 6 

Pfioton energy k/T 



QCD, 2-6 flavors 



QED: 1-2 species 

I , , I I I , , I , , , L 

0.2 0.4 0.6 0.£ 

Photon energy k/T 



FIG. 8. Left panel: Total constant under the log, equal to the sum of contributions Ctot{k/T) 
i log(2fc/r) + C2^2(fc/r) + Cbrem(fc/r) + Cannih(fc/r). This is the valuc which must be added to ln(r/moo) 



to obtain the total emission rate, as indicated in Eq. (6.1). The different curves (from bottom to top) are 
the same as in Fig. H. Right panel: Magnified view of the soft photon region. 
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Photon production rate Electromagnetic energy loss rate 




02468 10 02468 10 

Photon energy k/T Photon energy k/T 



FIG. 9. Total photon emission rate, together with the bremsstrahlung, inelastic pair annihilation and 
2 2 contributions, for two-flavor QCD with — 0.2. The left panel shows dT^/dk, divided by as Q^em T^, 
while the right panel shows rates weighted by photon energy. 

as (k/Ty^"^ but this asymptotic behavior turns out to set in only for k > 10 T. 

VII. CONCLUSIONS 

The above results represent the first complete leading order calculation of the photon 
emission rate from a hot, weakly coupled, thermalized quark gluon plasma at zero chemical 
potential. In addition to well known 2 2 particle processes, near-collinear bremsstrahlung 
and inelastic pair annihilation also make leading order contributions. A consistent treat- 
ment of these latter processes requires a rather complicated analysis in order to incorporate 
correctly the effects of multiple soft scatterings which may occur during the emission of the 
photon. These lead to an 0(1) suppression of these near-collinear processes. To account 
for this suppression, known as the LPM effect, we solved the integral equation, derived in 
Ref. [|TU|, which accounts for partial coherence between successive scattering events. 

For modest values of Ofg ~ 0.2, we find that near-collinear bremsstrahlung dominates 
the photo-emission rate for soft photons with energy k < 2T, and inelastic pair annihilation 
dominates the rate for very hard photons. A; ^ lOT. For intermediate values of k, which is 
the most easily detected range of photon energies, the 2^2 processes are of comparable 
importance to the sum of the near-collinear inelastic processes. The LPM suppression is 
substantial precisely where the inelastic mechanisms are dominant, but is fairly mild for 
intermediate values of photon energy. Parametrically, the soft photon emission rate behaves 
as dVy oc dk, which is less infrared singular, by a square root of k, from the result 
found in Refs. and [|T^ which neglect the LPM effect. 
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